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Abstract. By rewriting the Navier-Stokes equation in terms of differential 
forms we give a formulation which is abstracted and reproduced in a finite 
dimensional setting. We give two examples of these finite models and, in the 
latter case, prove some approximation results. Some useful properties of these 
finite models are derived. 



1. Introduction 

Increasingly the lines between algebraic topology, applied mathematics, differen- 
tial geometry, analysis, mathematical physics, etc., are becoming difficult to even 
imagine, as tools in one area become useful or essential in another. 

Tracing the history of one famous line of research, the study of the heat equation, 
one sees such an interplay. For example, the abstraction of the Euclidean analysis 
problem to the Riemannian setting led to spectral invariants and new notions of 
cohomology, while the further abstraction to an algebraic setting led to heat flow 
on graphs and simplicial complexes whose invariants include combinatorial counts 
such as spanning trees and their generalizations. At the same time, approximation 
theory was used to better understand mathematical invariants such as torsion, and 
conversely, techniques from algebra, topology, and analysis were (and still are) being 
used to construct models of heat flow. 

Fundamental in this story is some abstract formulation of the problem in ques- 
tion, so that tools in one area can be used to answer and pose questions in another. 

The purpose of this paper is to give one such abstraction for the Navier-Stokes 
equation, suitable for expression in both a smooth Riemannian manifold and flnite 
cochain complex setting. The results we obtain suggest a great potential for these 
two viewpoints to inform one another (e.g. Corollary |4]). Additionally, we indi- 
cate how this viewpoint may be used to simulate fluids, and expect it will lead to 
interesting dynamical and combinatorial problems in these flnite models. 

I would like to thank Jozef Dodziuk and Dennis Sullivan for helpful discussions 
related to this work. 



2. Viscous incompressible fluids 

In this section we explain how to express the Navier-Stokes equation in the 
language of differential forms on Riemannian manifold. The formulation obtained 
in the zero viscosity case, which corresponds to Euler's equation, agrees with W (p. 
588). We then give some basic properties which follow immediately from algebraic 
considerations. 
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The Navier-Stokes equation for a divergence free vector field representing an 
incompressible homogeneous fluid of viscosity > is given by 

(1) — + {u-V)u^ -Vp + uAu 

div u — 

where u is a time dependent vector field in R"^ and the pressure, p, is a time 
dependent function on M.^ . The latter equation, that the vector field u is divergence 
free, is referred to as the incompressibility condition. Using the vector identity 

{u ■ V)w ~ —u X curlu + ^V||m||^ 

we may re-write equation ([1]) as 

du 1 

(2) — = -li X curlw + -Vllwf - Vp + i/Au 

We now rewrite this equation in terms of a differential 1-form on M'^, thought 
of as a Riemannian manifold with the standard Euclidean metric, where 1-forms 
and vector fields are identified by the metric. Under this identification, the cross 
product is given by (w,r?) *(uj A 77), and the curl operator is ★d, where ★ is the 
Hodge-star operator. Also, the divergence free vector fields correspond to those 1- 
forms uj which satisfy d*uj = 0, and the Laplacian of a divergence free 1-form is given 
by ~d*dLd. Putting this together we can rewrite equation ([2]), the Navier-Stokes on 
the Riemannian manifold R^, as 

doj 1 

(3) — = - -k {oj A irdoj) + -dWuW^ - dp - lyd* doj d*u!^Q 

for a time dependent 1-form uj. For any time dependent 1-form uj(t) on a Riemann- 
ian manifold, we refer to this system as the Navier-Stokes equation in a Riemannian 
manifold. In the zero- viscosity case, u = 0, the equation is referred to as Euler's 
equation on a Riemannian manifold. 

Assuming our manifold is compactQ, we now make one more simplification. By 
the Hodge decomposition of differential forms we see that, since the left hand side 
of equation ([3]) is in the kernel of d* , the d-exact part of the right hand side must be 
zero. Therefore the Navier-Stokes equation in a Riemannian manifold is equivalent 
to 

duj 

(4) — ^ n -k (uj A -kduj)) - i^d* duj d*uj = 

where tt denotes the orthogonal projection onto the kernel of c?*, i.e. the co-closed 
1-forms. 

Setting J/ = we obtain the analogous Euler equation in a Riemannian manifold 

(5) ^ =7r(-*(wA*dw)) d*uj = 

From this formulation, several properties are apparent. First, if additionally ujq 
is d-closed (so uto is harmonic) then the right hand side of @ is zero, and therefore 
uj{t) = Wo is a steady state solution. Non-harmonic steady state solutions to Euler's 
equation are known to exist in general. 



^One may instead use differential forms with compact support, or those with rapid decay, as 
in the original Helmholtz decomposition. 
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Second, if uj{t) is a solution to (jl}, then by taking the inner product of both 
sides of Q with uj{t) we obtain 

(6) |||^f =.-2HM^f 

Here we have used that d* is the adjoint of d and that 

(7) (tt (— ★ (cj a *duj)) , cj) = (- ★ (w A -kduj),uj) — {diu, lu Aio) ~ 

In particular, we see that an Eulerian flow is norm preserving, while for a general 
Navier-Stokes flow the norm is non- increasing, as a function of the vorticity, duj. 



3. Algebraic and analytic considerations 

The equations of the previous section are well defined on the space of smooth 
differential forms. For the purposes of finite computation and approximation, one 
must work with more general spaces, e.g. piece-wise linear, piece-wise polynomial, 
or piece-wise smooth forms. As we now describe, there are particular difficulties 
with defining all of the operators d, A, and -k on these spaces or their generalizations. 

A particularly nice class of differential forms, sufficient for topology, are the 
Whitney flat forms |8], [6]. Roughly speaking, these are L°°-forms (bounded, mea- 
sureable, a.e.) whose exterior derivative is defined almost everywhere and is in 
A theorem of Wolfe, illuminated in Whitney's book [8], characterizes these forms as 
precisely the bounded linear functionals on chains with respect to a certain norm on 
chains. They are in fact defined for any Lipschitz manifold. These differential forms 
constitute a differential graded Banach algebra but, unfortunately, the Hodge star 
operator does not preserve this space. For example, the Hodge star operator does 
not preserve piece-wise linear, piece-wise polynomial, or piece-wise smooth forms 
whose singularities occur along the faces of a triangulation. 

On the other hand, the space of L^-forms does have a well defined and bounded 
Hodge star operator, and has the advantage of being a Hilbert space. We may even 
consider the space of i^-forms with d in L^, or the Sobolev space of forms with d 
and d* in . Unfortunately, the wedge product of forms is not well defined on any 
of these spaces. 

There is also the class of forms on an n-manifold given by fc-forms in L"/*"' whose 
exterior derivative is in ^^"/fe+i. By Holder's inequality, these form a differential 
graded Banach algebra, but again, the Hodge star operator is not well definecd. 

These apparent difficulties are addressed here by considering the space of 
forms as a module over the algebra of L°°-forms. Specifically, on a compact mani- 
fold, there is a continuous inclusion of Whitney fiat forms into the Hilbert space of 
L^-forms. If w is a Whitney flat form we can define T{uj) e by 

(8) {T{co), v) = (duJ, ^^V) for aU tj e 

where (, ) is the inner product on the Hilbert space of L2-forms. Indeed, w A77 G 
and dut G L°° C L^, and the right hand side defines for each uj a bounded linear 
functional of rj on (of norm at most a constant times ||dci;||2||a;||oo)- As we now 
explain, if uj happens to be a smooth form, this vector T{uj) is the time derivative 
of the Euler flow starting from uj. 



'It seems useful to have a conceptual explanation for these difficulties. 
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Let denote the bounded linear map from 1-forms given by left multiplica- 
tion by the 1-form lo G L°° . One can calculate that the adjoint of Li^, with respect 
to the inner product, is given by — L^i^. Therefore, we have 

(— * (oj A -kdu),!]) = {duj, uj Arj) for all rj £ L'^ 

From ([5]) and ([8]) we see that for a smooth time dependent 1-form w the projection 
of T(uj) onto co-closed 1-forms equals the time derivative of an Eulerian fluid 
flow, i.e. the Euler equation is given by 

(9) ^ = ^(Tiu:)) d*u; = 

where tt is the projection on the the kernel of d* . Similarly, we can re-express the 
Navier-Stokes equation, as 

(10) ^ = ^(T,H) d*oj = 
where T^(w) is defined by 

(11) (Ty (w), r/) = {duj, CO Arj) — v{dui, drf) 

The non-linear operator Ti, can be regarded as a section of the trivial bundle 
over Whitney flat forms whose fiber is L^-forms. The restriction of this section 
to smooth forms is an honest vector field, tangent to smooth forms, while for a 
general Whitney flat form it is not tangent. Of course, this makes solving the ODE 
starting from an initial Whitney flat form impossible. Nevertheless, this approach 
enlarges the space to one where the finite models for fluids of the next section can 
be formulated and compared by estimates. 



4. Finite models 

From the previous section we are motivated to consider the following structure: a 
finite dimensional cochain complex C" = {C^, 5} with graded commutative product 
U and positive definite inner product (, ) on C*. For any such structure, the inner 
product induces an isomorphism from to (C^)*, the dual of C^, and for any 
J/ > we can define a non-linear fiow on given by 

(12) |^ = vr(T.(c)) 

where tt is the projection onto the kernel of the adjoint 5* of (5, and T^{c) is defined 

by 

(r^(c), b) = ((5c, c U 6) - iy{Sc, Sb) 

Clearly, if S*Cq =0 then the flow starting from cq remains in the kernel of S* . 

A first example of this structure, associated to any simplicial complex, will be 
referred to as the toy-model. Let C* = C*{K) be the simplicial cochains of a 
closed simplicial complex K, with coboundary operator 6. There is a canonical 
basis of elementary cochains given by those whose value is one on a single simplex 
and zero elsewhere; in this way we can confuse an elementary cochain with the 
unique simplex on which it is supported. There is also a graded commutative non- 
associative product on C described easily in terms of the elementary cochains a, b 
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as follows: a U 5 is zero unless a and b intersect in exactly one vertex and span a 
(j + fc)-siniplex c, in which case, 

(13) aUb^eiaMjj^^c 

where e(a, b) is determined by 

orientation{a) ■ orientation{b) = e(a, b) ■ orientation{c) 

see [5] . Lastly, there is a positive definite inner product on C* defined by declaring 
the elementary cochains to be an orthonormal basis. 

A second example of this structure, which we call the Whitney model, is obtained 
from the triangulation of any closed Riemannian manifold. For this we take C* to 
be simplicial cochains of the triangulation, U to be the same product as in the last 
example, and (, ) to be the Whitney metric which we now describe. 

There is a injective cochain map W from C* to the space of Whitney flat forms 
[8] . Indeed, W is defined on an elementary cochain a supported on a single simplex 
a = [pojPIj • ■ • tPj] with barycentric coordinates corresponding to the i^^ vertex 
Pi of a by 





W{a) = jl d^Q A • • • A d/ij A • • • A dfij. 

i=Q 

and then is extended linearly to all of C* . Moreover, RW — Id where R is the 
integration map from forms to cochains of K . We will not indicate the role of K 
in the maps W or R. 

We define an inner product on C*, as in |3^, by the restriction of the inner 
product to the image of W, and denote it by the same notation: (a, 6) = {Wa, Wb). 
It is shown in [3] that this is a positive definite inner product. 

The ingredients of this Whitney model have nice approximation properties, as 
indicated in the following theorem, proved in P|. 

Theorem 1. Let K be a triangulation of a manifold M with mesh rj. There exist 
a constant C and positive integer m, independent of K such that 

\\W{R(jJi U RUJ2) - cji A W2II < C • A(wi,a;2) • 11 

where 

A(c^i,c^2) = llc^illoo • Arc^2|| + ||c^2||oo ' + Arc^i|| 

for all smooth forms loi,lo2 G 0,{M), where \\ \\, \\ \\oo are the and norms, 
respectively. 

Remark 2. Applying the previous theorem to the constant function a;2 = 1 shows 
that \\WRuj — uj\\ < Crj, a result previously obtained in [3]. 

This theorem makes precise the sense in which, in a sequence of triangulations 
with mesh converging to zero, the cup product converges to the wedge product. 
As is always the case, in a sequence of refining triangulations we require that the 
shapes do not become too thin, i.e. the fullness must be bounded; see for example 
[5] for details. We call this a nice sequence of subdivisions. 

In light of these convergence results, one might hope to prove a strong approx- 
imation property of the Whitney model, like the following. For a solution ijj{t) to 
the Navier-Stokes equations, we could ask that for any given degree of closeness 
e > 0, and desired time > 0, one can choose a small enough mesh so that the 
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solution c(t) — (i?ajo)(t) to the Whitney model starting from Rujo is within e of 
uj{t) for alio < t < S. This is in fact too strong of a request since, for example, the 
equation may not be stable near a given smooth initial condition. 

Instead we are able to show that the local operators used to define the Navier- 
Stokes are approximated well, in a weak sense, by the Whitney model. In particular, 
the Whitney model provides a weak approximation to the time derivative of the 
Navier Stokes flow starting from a smooth solution. 

Theorem 3. Let lo be a smooth form on a closed Riemannian manifold M and 
Ki,K2, ... be a nice sequence of triangulations of M . Let c = Ruj for each trian- 
gulation. 

Let IT {T^Lo) and tt {T,jc) be the tangent vectors defining the Navier-Stokes flow on 
M and the Whitney models, respectively, defined as above. Then Wtt(Ti,c) converges 
weakly to 7r(T^w) in the sense that, for all rj G L^ , 

as the mesh tends to zero. 

Proof. We consider the two summands of separately. For all cj and rj we have 
that 

{WSRuj,rj) = {WRdLo,ri) {duj,rj) 

by Remark [21 For the non- linear term we denote Tq by T, and have that 

\im{WTTTRu, Tj) ^ lim(I^7rri?cj, WRrj) = Y\uy{TRuj, ttRt]) 

by Remark[2]and the fact that orthogonal projection is self adjoint. It is known that 
the Hodge decomposition of cochains, with respect to the Whitney metric, converges 
to the smooth Hodge decomposition (see [1] Theorem 2.10, and [2] Lemma 3.18 for 
a recent improved estimate). Therefore, we have that linnrRr] = HmRirr], and that 
the previous expression is equal to 

lim(Ti?a;, Rttj]) = \im{dRuj, Ruj U Rttt]) = \im{W5RLu, W {Ruj U Rttt])) 

Now we use that the cup product converges to the wedge product, Theorem [TJ so 
the previous expression is equal to 

(dw, w A TTjy) = {Tuj,Trri) = {-nTuj^Tf), 

and this completes the proof. 

We remark that the rate of convergence is on the order of mesh size constant 
factor depending only on the norm of uj, rj, their derivatives, and some universal 
constants independent of the triangulations. □ 

There is a purely algebraic viewpoint of the argument above. Note that the 
operator -kL^^-k corresponds d*-homologically to the cap product of uj on the complex 
of forms (with differential d*). This defines a differential module of {yi,d*) over 
(r2(M), d, a). By the argument above, the analogous capping operation on C* , 
defined as the adjoint of cup product, converges weakly to this module structure 
on forms. One can of course state the flow purely in terms of this cap product, 
though we have not chosen to do so. 

From the previous result we also obtain the following result concerning steady 
state solutions. 
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Corollary 4. Let M he a closed Riemannian manifold. The following is a necessary 
condition for a smooth 1-form uj to be a steady state solution to the Navier-Stokes 
equation: for every nice sequence of subdivisions of a given triangulation of M the 
sequence of cochains c = Rlo given by integrating uj over each subdivision must con- 
verge weakly to a steady state solution of the Whitney model as the mesh converges 
to zero. 

Proof. Steady state solutions are those forms w which satisfy tt {T^{ijj)) = 0. By 
the previous proposition, this is approximated weakly by Wt:{T^c). □ 

We remark that for the Whitney model associated to a particular triangulation, 
the integral of a smooth steady state solution may not be a steady state solution in 
the Whitney model. Also, at any finite stage, the Whitney model may have steady 
state solutions that are not given by the integral of smooth steady state solutions. 
For example, it is clear from (|12l) that a cochain c in the Whitney model is a steady 
state solution to Euler's equation if the "vorticity" 6c £ is orthogonal to the 
image of the operator Lc : ^ given by cupping with c. This condition 
is intimately related to not only the metric, but also to the combinatorics of the 
triangulation. 

So far we have considered only piece-wise linear differential forms for our finite 
models. It is clear that using higher degree polynomial forms would yield better 
approximation results. For example, the recent finite element spaces in [2] may be 
used in a similar way. 



The finite models of the previous section, consisting of a finite cochain complex 
C* — {C^ ,6} with graded commutative product U and positive definite inner prod- 
uct (, ), also enjoy the properties derived in section 1. In particular, the harmonic 
elements, i.e. those satisfying 6c — 6*c = 0, are steady-state solutions to ([T2l) . Also, 
the same proof as in ([6]), now using the graded commutativity of U, shows that 



holds for any co-closed solution to (|12p . 

Since the finite models are defined by algebraic equations on a finite dimensional 
vector space, the flow for (fT2|) is defined for all time, until blow up. From ()14|) we 
see that any finite model Navier-Stokes or Euler flow is in fact defined for all time 
since the norm is non-increasing. 

We remark that the toy-model and Whitney model may be implemented com- 
putationally by expressing the coboundary operator 6, the cup product U and the 
inner product (, ) in terms of the basis of elementary cochains. There is a local 
calculation that one must perform to determine the Whitney metric from the given 
Riemannian metric. The last required operator is the orthogonal projection onto 
the kernel of 5* , which may be computed for example by first computing the Hodge 
decomposition. 

It is worth noting that the Whitney metric and cup product (|13l) are semi-local 
with respect to the basis of elementary cochains: two elementary cochains have 
non-zero (inner) product only if the they are supported on simplices which are 
faces of a common top dimensional simplex. The inner product of the toy-model 



5. Remarks on finite models 



(14) 




2u\\dc\\^ 
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is even more local, as the elementary cochains are an orthonormal basis. On the 
other hand, the orthogonal projection operator in ([T2l) is non-local. 

We do not presently have a means to compare the two finite models we've given, 
though it is conceivable that the toy model alone will produce realistic fluid models. 
To make this precise, one needs a notion of a (potentially soft) map or morphisni 
between flnite models or sequences of them. In [7 a related issue is addressed for 
the algebraic structure of Euler's equation. We hope the pursuit of similar algebraic 
ideas along with the formulations of this paper will lead to a further understanding 
of fluids and fluid models. 
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